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A two-dimensional Fourier transform of hadron form factors allows to determine their charge 
density in transverse space. We show that this method can be applied to any virtual photon induced 
transition, such as y*(q) + N — > nN. Only Fock states that are common to the initial and final states 
contribute to the amplitudes, which are determined by the overlap of the corresponding light-front 
wave functions. Their transverse extent may be studied as a function of the final state configuration, 
allowing qualitatively new insight into strong interaction dynamics. Fourier transforming the cross 
section (rather than the amplitude) gives the distribution of the transverse distance between the 
virtual photon interaction vertices in the scattering amplitude and its complex conjugate. While 
the measurement of parton distributions in longitudinal momentum depends on the leading twist 
approximation (— q 2 — > oo limit), all q 2 < values contribute to the Fourier transform, with the 
transverse resolution increasing with the available range in q 2 . We illustrate the method using QED 
amplitudes. 

PACS numbers: 13.60.-r 



I. INTRODUCTION 

The photon has a pointlike coupling to quarks, which makes it a valuable probe of hadronic processes at any 
virtuality q 2 . The nucleon charge radius is determined from the slope of eN — > eN as q — > [1]. In this limit 
the nucleon acts as a static source. At higher q 2 one expects to be able to map the charge distribution in more 
detail. However, the quarks are highly relativistic and move as fast as the probe. Hence the time difference between 
spatially separated photon interactions cannot be neglected. The standard three-dimensional Fourier transform, 
which is appropriate for non-relativistic systems such as atoms, assumes an instantaneous photon interaction and is 
not applicable for tracing quarks in hadrons. 

The precise relation between the spatial distribution of quarks and hadron form factors was long obscure. It was 
uncovered only recently via a somewhat circuitous route involving deep inelastic scattering, eA^ — > eX (DIS) [2- 
6]. In the Bjorken limit (q 2 — > — oo) the DIS cross section is dominated by photon scattering off the same quark 
in the amplitude and its complex conjugate, with the interaction vertices separated by a light-like distance. The 
distribution of the quark longitudinal momentum fraction x is therefore most concisely expressed in terms of nucleon 
wave functions defined at equal Light-Front (LF) time x + = t + z, [2, 7, 8] 
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Here ip~* (xi, hi) is the wave function of a nucleon Fock state with n quarks and gluons, all at the same x + and with 
each constituent carrying momentum fraction Xi and transverse momentum ki relative to the nucleon. 1 A Fock state 
only contributes to fq/i\r(x, Q 2 ) provided it contains a quark with Xk = x. The absolute square of the wave function 
allows to regard the quark distribution as a probability density. 2 

The possibility to access generalized parton distributions (GPDs) in deeply virtual Compton scattering, eN — > ejN' , 
rekindled interest in studying the distribution of quarks also in transverse space (impact parameter) [4-6]. Fourier 
transforming the transverse momentum p(N')—p(N) dependence of the GPD was shown to give the impact parameter 
b distribution of quarks (after an extrapolation to p + (N) = p + (N')). The extraction of GPD's from scattering data is 
demanding, hence our knowledge of quark distributions q(x, b) in both longitudinal momentum and impact parameter 
remains model dependent. 



1 Throughout this paper we use bold symbols to indicate 2-dimensional transverse vectors. Inessential helicities are suppressed. 

2 This is only approximate [9], since Coulomb scattering of the quark between the photon vertices is neglected in (1). 
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The GPD's reduce to electromagnetic form factors when integrated over the longitudinal momentum fraction, 
schematically J dx q(x,b) — F(b). This revealed the desired relation between hadron form factors in momentum 
space and their charge density in transverse space. Form factors are much easier to measure than GPD's, which 
allowed to plot the nucleon density distributions without model dependence [10, 11]. In form factors the momentum 
transfer to the target equals the photon momentum, q — p(N') ~ p(N). There is no notion of a "leading twist" 
approximation for form factors. The Fourier transform wrt. q gives the charge density of quarks as a function of 
impact parameter b. For the Dirac form factor (Q 2 ) of the nucleon, 3 

po(b) ee J^e-^ b ^(N(p+,y)\J+(0)\N(p+,-lq))=J™^QJ (bQ)F 1 N (Q 2 ) 
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Here the (xi,bi) are LF wave functions of a nucleon state \N(p + ,b^ = 0)) with 'plus' momentum p + = p° + p 3 
and transverse 'center-of-momentum' b^ — Y^ x i^i = 0. The quarks and gluons in each n-parton Fock state have 
longitudinal momenta — xip + and impact parameters 6^. Only quarks at transverse position b^ = b contribute 
to the charge density at b. The quark distribution in impact parameter (2) is analogous and complementary to the 
standard parton distribution f q /N(x,Q 2 ) (1) in longitudinal momentum. 4 

An expansion similar to (2) pertains also for transition form factors measured in eiV — > eN* . As we recall below 
(Section II) the expression remains diagonal in the LF Fock basis: only Fock states that are common to N and N* 
contribute. The corresponding wave functions ij;^ , tpn being distinct their product is no longer positive definite. 
Nevertheless, the impact parameter distribution reflects the transverse size of the transition process, and has been 
studied using data on several nucleon resonances [11, 12]. 

The expression (2) for the impact parameter distribution only assumes the general LF Fock expansion of the intial 
and final hadronic states. Hence it can be applied also to states with several hadrons in the final (and initial) 
state. This allows to study the transverse size of photon scattering processes as a function of the relative momenta 
of the final state hadrons. The method can even be applied at the level of cross-sections, thus not requiring a 
knowledge of the phase of the scattering amplitudes in the Fourier transform. One obtains then the distribution of 
the transverse distance between the photon vertices in the amplitude and its complex conjugate. Since no leading 
twist approximation is implied this type of analysis is particularly suitable for data at moderate values of q 2 , provided 
only that the contribution of the J + current can be isolated. The resolution in impact parameter improves with the 
range of q 2 for which data is available. 



II. BASIC FORMALISM 



We first recall [4, 5] the basic steps which lead to the expression corresponding to the nucleon density (2) for any 
final state /. The lepton scattering amplitude is 

M(IN -> t'f) = -e 2 u(£')^u(£)^ J d 4 xe~^(f\J»{x)\N(p)) (3) 

where q — £ — £' is the virtual photon momentum. Using LF spinors [7] quantized along the negative z-axis and 
neglecting the lepton mass, 

u(l', X e )^u{£, Xg) = . * [2tTlp - £-q„ - q-lp + n^£ ■ q + zA.e^n^V] (4) 

y , £-{£--q-) 

where the light-like vector n = (2 + , 0~, OjJ satisfies n ■ p = p~ = p° — p 3 for any vector p. The J + hadron matrix 
element dominates in (3) in the high energy limit, £~ — > oo at fixed q. This limit is also required to formally access 
all momentum transfers q. Hence we consider 

M(eN^£ / f)^-e 2e -(f(p f )\J+(0)\N(p))(27r) 4 S i (p f -p- q ) (f ^ oo) (5) 



3 We follow the conventions of Ref. [5]. 

4 The form factor has a single photon vertex and hence no Wilson line. The relation (2) is exact (up to higher order electromagnetic 
corrections) insofar as the LF Fock expansion of hadrons is exact (contributions from partons with Xj = are neglected). 
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The matrix element (f(pf)\ J + (0)\N(p)) may be viewed as a generalized form factor. Apart from its momentum 
Pf = p + q there is no restriction on the final hadronic state / which could, e.g., consist of many hadrons (see Section 
III). 

The J + quark current projects on the q+ component of the quark field, 



J + (x) = e q q(x)^ + q(x) = 2e q q\_ (x)q + (x) 



(6) 



where q+(x) 



:) and the light-like vector n = (0 + , 2 , 0j_) satisfies n-p = p + for an arbitrary 4- vector p. The 



q+ quark field may be expanded in LF creation and destruction operators at a given LF time. For x + = 

dk+d 2 k 



q+(0 + ,a; ,x) 



16tt 3 /c+ 



9(k+) b(k+,k)u + (k+)e~ l i k x ~ +tk ' x + d\k+, k)v + (k+)e 



(7) 



The u+ and w+ spinors are independent of the transverse momentum k and normalized according to u\_(k + )u + (k + ) = 
k + . In terms of operators at a fixed transverse position x, 



b(k+,x) = 



d 2 k 
16tt 3 



3 lk - x b{k + ,k) 



the quark field is expressed more simply as 



ydk + 
— 8(k + )\b(k + ,x)u + (k+)e- l 2 kx ~ + d 1f (k + ,x)v + (k+)c 

The transverse momentum eigenstates may be expanded in impact parameter states 

\p+,p) =4tt J d 2 be lp b \p+,b) 

which have the LF (x + = 0) Fock expansion 



(8) 



(9) 
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inn 
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The n operators in each Fock state create quarks (tf), antiquarks (cP ) and gluons (a*) with longitudinal momenta XiP + 
at transverse positions ?v The specific advantage of the LF Fock expansion is that a hadron with any longitudinal 
momentum p + and transverse position b is described by the same LF wave functions ij) n (xi, bi — 6), which depend 
only on the relative coordinates of the partons. 

The quark field (9) eliminates an operator b^(xkP + , bfe) at b^ = x from the Fock expansion (11), according to the 
anti-commutation relation 



\b(k + , b),bUk' + , b')\ = ^-k+6(k+ - k' + ) S 2 (b - b') 



(12) 



Thus, suppressing the contribution of the creation operator d'(k + ,x) (see below) 

_ > ,rr r 1 dx 

(47T) 



1 

q+ (0+,x-,x)\p+,b) = T^Efn 



4nd 2 bi 



(5(1 - x i) §2 i b - X! Xlb ^ ^n(x t , bi - b) 



Y J [{-^ Pk 5\ b k-x)u + {x W + )e- 1 ^ + ^\[b\x l p + M)d\ W( )] |0) 

k i^k 



(13) 



where the sign (— \) Pk related to operator ordering will be irrelevant, since according to (6) the J + matrix element in 
(5) is the overlap of two states of the form (13). 

To allow a simple interpretation of the amplitude (5) it is essential to choose a frame where p^ = p + . 5 A photon 
with q + = cannot create a qq pair, causing the matrix element to be diagonal in the number of incoming and 



5 In the case of GPD's this condition implies an extrapolation from the experimentally accessible kinematic region. For form factors it 
amounts to a choice of frame. 
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outgoing quarks. In fact, the initial and final Fock states are identical. As seen from (13) the J + (0) current interacts 
with a single quark or antiquark 6 at = 0± in \N), and similarly in (/|. The remaining n— 1 partons in \N) must 
thus be identical to those in (/|. The constraints J2i x i = 1 m the initial and final states forces also the momentum 
fraction Xk of the struck quark to be the same. The "center of momentum" constraint b = J^ajjbj in (13) requires 
the impact parameters of the initial and final states to be equal, 

^(f(p + ,b f )\J+(0)\N(p + ,b N )) ee *p(bf-b N )A fN (-b N ) (14) 
where, after a shift of integration variables bi — > bi + bjv, 

AfN(b) = 7^[f[ / dx * /47rd 2 6 l ]5(l-^x05 2 (^x l 6 4 )'/'r(^,^)^(^,^)5Ie fc «5 2 (b fe -b) (15) 

n i=l JO J i i k 

This expression for the current matrix element in impact parameter space is central for the applications we consider 
below. For / = N the positivity of \tp^ (xi, bj)| 2 allows the Fourier transform (2) of the elastic (helicity non-flip) form 
factor to be interpreted as a charge density. Even when the final state differs from the initial one its electro-excitation 
still proceeds only via Fock components n which are common to both. 

As already indicated in (2), the Fourier transform wrt. q of the generalized form factor in (5) should be done in a 
frame where the nucleon and photon momenta are 

P = (p + ,P~,-kl) 

q = (0 + ,q-,q) (16) 

Pf = (p + ,p~ + \q) 



The excitation amplitude in impact parameter space is then, using (10) and (14), 

: {f(p f )\J+(0)\N(p)) = (17) 



rf2 9 r -ic b 1 
(2tt) 2 2p+ 

= 1 10f d2bNd2b f ^• (b+ ^ N+ ^ /)1 0r(/(p + ^/)l^ + (o)l^b + ^ Ar ))=^ /A ,(b) 

The expansion (15) shows that A/N(b) gets contributions from LF Fock states that are common to the initial and 
final states (localized at bjy = 6/ = 0) which have a quark or antiquark at transverse position b/, = b. The range of 
AfN(b) in b thus reflects the transverse size of the transition process. 

The above analysis has previously been applied to elastic and transition electromagnetic form factors [10-12]. The 
Fock expansion (11) is, however, completely general and applies also to states |/) that consist of several hadrons. This 
makes it possible to measure the transverse shape of the hadronic states that contribute to 7* + i — >• / transitions, 
for any states i and /. 



III. TWO-BODY FINAL STATES 



The momentum pj = p + q of the final state / varies with q in the Fourier transform (17), hence the p/-dependence 
of all Fock state wave functions in the expansion of \f(pf)) must be known. As seen from (11) the LF wave functions 
depend only on the relative coordinates of the constituents, not on the total momentum of the state. Final states 
I/) = \h%, . . . , h n ) consisting of several hadrons may be regarded as a particular type of hadronic state, where we are 
free to specify the relative momenta of the hadrons, each one of which has its own (non-perturbative) Fock expansion. 
The multi-hadron Fock amplitudes must conform with the general LF rules to ensure the frame independence of the 
state I/). In this Section we specify the LF Fock expansion and the Fourier transform for a two-body (irN) state, and 
illustrate it with a tree-level QED amplitude. The multi-hadron case is considered in Section IV, where we discuss 
the Fourier transform of the cross section. 



Due to the anti-commutation of the d-operators the charge in (15) has opposite sign for quarks and antiquarks. 
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A. Transverse shape analysis of ~y* N — »■ ttN 



The standard LF Fock expansion in transverse momentum space for a single pion is [5, 7] 

\n(p+, Pl )) = 167r 3 ]r[n f ^]s(l~Y,^) S \T, k ^n(^k l )f[b\x l pt,x lPl +k l )---\0) (18) 

where • • • stands for the operators which create the remaining n — 1 partons of the Fock state. As noted above, the 
wave functions ip^Xi^ki) are independent of the pion momentum p±. The 'phis' momentum of parton i is Xip\ and 
its transverse momentum is XiP 1 + ki. The restrictions on the Xi and ki implied by (18) ensure that the parton 
momenta sum to the total pion momentum in each Fock state. 
For a 7riV state we have then the double expansion 



(19) 



x 5- 



I'' TV UN 
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which should be transformed into the standard form (18), where parton momenta refer to the total momentum 
Pf = pi + p2 of the state. We parametrize the pion and nucleon momenta in terms of a momentum fraction x and 
relative transverse momentum k, 



Pi = XPf 

P + = (i - x)pi 



p 2 = (1 - x)p f - k 



(20) 



(21) 



where p~j = p + and p^ = |q in the frame (16). The momentum fractions of the pion and nucleon constituents wrt. 
p + are then x\ = xxi and y'j = (1 — x)yj, respectively. Using this and integrating over x gives 

/ dxs{i - J2 xi)s(i - vi) = »(i - - E < E ^ 

i j i 

where x = ^\ x\ on the rhs. The transverse momenta of the partons may be expressed as 

XiP^ ~\~ ki = x^p y ~f~ Aj,^ Ai^ Aj^ ~f~ k x ^ I x 

.'/,/>, +Zj = V'jPf +4 *j = f-i ;! - X) 

which gives 

% j i j 

For an |/) = \tcN) state specified in terms of the relative hadron momentum by a wave function ty* {x, k) we get 



(22) 



(23) 



\*N(p+,p f ;*f)) 



dx 



d 2 k 



o y/x(l-x)J 16tt 3 

1 dx'i dy'j f d 2 k\ d 2 t' 3 



16tt 



e [n 



&(x,k)\7r(p 1 )N(p 2 ))= (24) 

<i-E*i-E^> 2 (E fe/ *+E4 



167T 3 167T 3 



X X[l — 



(i - x)*/(x,fc)^(^ f K - x ^t{^ 4 + Y^- X k ) 

II [^(x^,x' iPf + *o • • •] [n , =&w> + 4) 



(25) 



() 



where x = J^i x '% an d ^ = J2i on the rhs. This Fock expansion has the standard LF form, implying that the 
superposition of ttN plane wave states should be given by (#, k) with x and fc determined by the relations (20) in 
any frame. 7 The standard normalization condition (suppressing the helicities) 



(irN(p' + ,p'^ f )\TTN{p + ,p;^ f )} = 16tt 3 p + 6(p + - p' + )6 2 (p - p') 



implies 



dx 



d 2 k 
16tt 3 



\^ f {x,k)f 



(26) 



(27) 



For the wave function fyf(x, k) to preserve the invariant mass of the irN state it should have support only at fixed 



p) = (Pi +P2) 2 



x 



1 — x x{l — x) 



(28) 



This would be satisfied, e.g., by standard partial wave analyses. A superposition defined by the appropriate spherical 
harmonics in the rest frame {p^ = 0) determines directly the frame independent wave function ^>^(x, k). However, 
it is not necessary to constrain the mass (28) to be fixed. States with different mass that are produced at the same 
q + = and q will differ wrt. q~ , which does not affect the Fourier transform. 

It is instructive to express the ttN states also in impact parameter space, again following the conventions for partonic 
states. For the state (24) 



\nN(p + ,p;^ f )) = 4tt J d 2 b e lp ' b \iTN(p + , b ; * / )) 



(29) 



where 

\nN(p + ,b- 1 ^ f )) = 4tt 



dx d 2 bn d 2 bM 2 



\J x(l — x 

The hadronic wave functions are related according to 



S 2 [xb n + (1 - x)b N ] ¥ (x, 6 W ) \ir(x P +,K + b)N((l - x)p+,b N + b)) (30) 



d 2 k 
16tt 3 



exp 



ik ■ b n 

1 - x 



ty f (x, k) 



The normalization condition in impact parameter space is 

47r J dx d 2 K d 2 b N S 2 [xK + (1 - x)b N ] \^(x, b„) 



(31) 



(32) 



B. Illustration: IN I' N* (1440) I'ttN 

We illustrate our approach by considering the Pll(1440) (Roper) resonance contribution to the 7*iV — > irN ampli- 
tude. The charge density given by the TV — > iV*(1440) transition form factors was previously determined in Ref. [12]. 
Here we include the resonance decay, N* —¥ ttN, with a fixed relative momentum between the pion and the nucleon. 
Thus the hadronic wave function {x, k) in (24) is a (5-function in x and k, compatible with the constraint (28). 

Connecting the N ^> N* matrix element given in Eq. (1) of [12] with the N* — > irN decay amplitude using the N* 
propagator we get the desired contribution to (n(pi)N(p 2 , A 2 )| J + (0)\N(p, X))/2p + at p 2 — M 2 , 



2p~ 



J2«Pl)N(p2,\2)\N*(p f ,\f)) 



M 2 



M 



— ^2,^)75—^ 



2 P 



iMT 



(N*(p f ,\ f )\J + (0)\N(p,\)) 



NN" 



M + m N 



u{p, A) 



(33) 



7 The Fock expansion (25) is for the non-interacting \irN)out state at x + — > 00. The tpn wave functions in (15) describe the state (24) 
at x + = 0, before the pion and nucleon have formed. 
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We evaluate this amplitude in the frame (16) where Q 2 = —q 2 = q 2 . The final state momenta p\ and P2 are 
parametrized as in (20) and we use the LF helicity spinors given in [7]. The N* and nucleon masses are denoted M 
and tojvj respectively, T is the total width of the N* and g* is the irNN* coupling constant. The result for the spin 
flip and non-flip amplitudes are 



A^L^k) = r- — — \ k^Fr'(Q 2 ) + Qe^ M{ \. l ] mN F? N \Q 2 ) j> (84) 
Vl — xMT { M + m N J 

where the transverse momenta are expressed as q = g_i_(cos0 g ,sin</> g ) and k — fcj_(cos<^.,sin^>fe). 
Fourier transforming to impact parameter space as in (17) gives 

A\ N + {b;x,k) = VY ^ MT { [M(l -x)- m N ]p$ N \b) + ik^*" N * (b)} 

Af_(b;x,k) = VT ^ MT {k^p% N \b) -i[M(l -x)- m N ]e^pf N \b)} (35) 
where b = 6(cos ^>j,, sin </>;,) and the charge density distributions 

NN* fi\ f°° T t^NN" f/~i2\ 



p NN* {b)= rdQQ^J l(bQ)F NN- m m 

J 2tt M + m N 

are expressed in terms of the transition form factors similarly as in Ref. [12]. The impact parameter amplitudes (35) 
have the Fock expansion given in (15). 

When one does not assume a specific intermediate irN state the 7*iV — > irN amplitude can be expressed in terms of 
six invariant amplitudes specified, e.g., in [13]. The LF helicity amplitudes in the frame (16) can be evaluated in terms 
of such invariant amplitudes similarly as above, allowing to determine charge densities using existing parametrisations 
of experimental data. 



C. Illustration: Ifi £' fif 

We denote the photon matrix element in the amplitude (5) for l/i — ¥ t' jjrf by 



^7,a 2 = ^ T (/i(Pi,A 1 ) 7 (P2,A 2 )|J + (0)| A1 b,A = 1)) 



(37) 



where p + q = pi + P2 — p / and (as indicated) the initial muon has helicity A = \ . At lowest order, using LF helicity 
spinors [7] in the frame (16) and the parametrization (20), diagrams (a) and (b) of Fig. 1 give, respectively, 

e • k e_ • [k — (1 — x)q] 



(1 - x) 2 m 2 +k 2 (1 - x) 2 m 2 + [k - (1 - x)qf 



(38) 



where e\ ■ k = — \e lX ^ k |fe|/v2- The corresponding expressions for the other helicity amplitudes are given in the 
Appendix. The Fourier transform (17) gives the amplitude for the virtual photon to interact with a muon at impact 
parameter b, when the center-of-momentum of the initial and final states is at zero impact parameter: 



A^i + l{b-,x,k)=2e^fx~ 



■ k 



(1 — x) 2 m 2 + k 



;S 2 (b)~ 



-i<t>b 



2V2n 1 



exp 



kb 

1-x 



(39) 



where b = b(cos sin fa). The first term in (38) arises from the diagram of Fig. 1(a), where the virtual photon 
vertex is before the real photon vertex on the muon line, i.e., the exchanged photon interacts with the initial muon. 
As expected, it contributes to (39) at the initial impact parameter b = 0. The second term is generated by Fig. 1(b), 
where the virtual photon interacts with the muon after the emission of the real photon, and the 6-dependence reflects 
the impact parameter distribution of the final state muon. 




(a) 




(b) 



FIG. 1: The two Feynman diagrams contributing to the QED process £fi — ¥ i'llj, when the photon is emitted from the muon. 



The amplitude (39) is in fact given precisely by the overlap (15) of the LF Fock amplitudes of the initial (/i) and 
final (fij) states. According to (37) the wave function ^(x, k) which describes the final state as in (24) is a <5-function 
in x and k. In impact parameter space (31) gives 



*(x', b) = S(x' - x)Jx{\ - x) exp (i ^ & 

\ 1 — x 



(40) 

The first term in (39) corresponds in (15) to the single particle (n = 1, initial muon) Fock state contribution, which has 
Xn = 1 and 6^ = 0. The coefficient of S 2 (b) in (39) must therefore be the c.c. of the single muon wave function in the 
fj/y final state, "0*(/i7 — > fi), which is the same as ip(p — > /ij) (with reversed sign due to the LF energy denominator). 
Multiplying the /i — > fij LF wave function given in, e.g., [14], 



2e 



il>\ x+1 (x,k) = 

+ y/ 1 — X (1 — X) 



k 



k 1 



(41) 



by ^*(x' , b) and integrating over x' indeed reproduces the coefficient of 5 2 {b) in the first term of (39). The second 
term in (39) arises from the n = 2 (/17) Fock state contribution in (15). It is readily seen to be the product of 
47r/(l — x) 2 arising from the integrations in (15), the final state amplitude (x', b) and the [i — > ^7 LF wave function 
in impact parameter space [15], 



4\/27r 



-emyj\ — xe % ^ b K\{mb) 



(42) 



Given the explicit expression for the QED amplitude (37) we may also consider ^7 final states with fixed impact 
parameter b'^ of the final muon. Choosing 

V(x',b') = 8(x'-x)y/x(l-x) ^^- 5 2 (b' -b') (43) 

47T 

we find from (31) 

*(x',fc) = 5{x - x)Jx(l- x) exp ( - i^^) (44) 

V 1 — x J 

Integrating over the relative momentum k of the final state with weight $>*(x', k) according to (24), 



d 2 k k- b'„ 

exp 1- 



16tt 3 



^ — em\fx(\ — x) e % ^ h '^ K\ tmb' 

AV2tt 2 V ; V ■ 



-1 



The Fourier transform (17) in q then gives 



-S^(b) + 6^(b-b') 



(45) 



(46) 



Thus the virtual photon interacts either with the initial muon at b = or the final muon at b = b'^. In accordance 
with (15) the distribution is determined by the LF wave function (42) for /1 — > fij (with the sign change noted above). 
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IV. CROSS SECTIONS IN IMPACT PARAMETER SPACE 



The superposition (24) and Fourier transform (17) discussed above require a knowledge of the phase of the scattering 
amplitude (f(pf)\ J + (0)\N(p)). Since a partial wave analysis is practical only for a limited subset of all amplitudes it 
is interesting to ask what information about the transverse structure of the scattering process can be obtained from 
a Fourier transform of the measured cross section. As we next discuss, this gives the distribution of the transverse 
distance between the photon interaction vertices in the amplitude and its complex conjugate. 

As in the case of the amplitude (5) we need to isolate the contribution of the J + current. Here we again consider 
the high energy limit s ~ £~p + — > oo at fixed momentum transfer q — £ — £'. The Lorentz invariant cross section can 
then be expressed as 



. da(W -> ej) 
dq~ d 2 q 



7r q 1 



dU 



f 



+ (f(p f )\J + (o)Wp)} 



(47) 



where dUf is the phase space element of the hadrons in /. The frame (16) can be reached from the IN CM by 
a rotation SO ~ |q|/^ _ around the normal to the lepton scattering plane. In the £~ — > oo limit the rotation is 
infinitesimal and does not affect the finite momentum transfer q. Then the Fourier transformation below can be done 
directly in the £N CM. 

For a state / with N^ hadrons of momenta pi , 



dU f (N h ) = 



N h 

n 



dpj d z p, 
\ (2^) 3 2^ 



With a LF parametrization as in (20), 



XiP~f 



Pi = %iPf 



(2^) 4 <5 4 (p + g-5>) 



,N h ) 



where Pf = X^jPii we obtain 



<m f (N h ) = 



2(2tt) 4 



'N h 

n 



dX'i d ki 



■Pf) 



(48) 



(49) 



(50) 



The initial nucleon N and final state / in the matrix element of (47) may be Fourier transformed (10) in the frame 
(16), where p f = p = \q and q + = 0. According to (14) the matrix element is diagonal in impact parameter. Thus 



S fN (b) 



d2 1 e -<«-i 
(2tt)2 



2p 



: (f(p f )\J + (0)\N(p)) 



= fd* bq AMb q )A} N ( bq - b ) 



(51) 



Altogether we get for the Fourier transformed cross section, 



d 2 q 

(2tt)2 



*f) 



d 2 q 



f 



n 



dx>i d k<i 
(2^)322;, 



(52) 



As indicated, the cross section may include several final states / with different hadron multiplicities Nh- The am- 
plitudes AfN{b q ) defined by (17) can according to (15) be expanded in terms of Fock states common to N and /. 
With the initial and final states located at zero impact parameter the struck quark is at impact parameter b q . Hence 
Sfjq{b) gives the distribution in transverse distance b between the quark struck in the amplitude and in its complex 
conjugate. It has a real part that is even under b — > b and an imaginary part that is odd. A non-vanishing imaginary 
part requires that the squared matrix element in (51) changes when q — > —q. This can be caused by a correlation 
between q and a transverse direction defined by the final state /. For example, in (38) the amplitude depends on the 
angle between q and the relative transverse momentum k of the muon and the photon in the final state. A direction 
can also be specified by a transverse polarization in the initial or final state. 

The final phase space integral in (52) refers to the internal momenta of the final state /. E.g., in the particular 
case of |/) = \ir(pi)N(p2)) , with p\ and P2 defined by (20) and the hadronic wave function (cc, /c) chosen to be a 
(^-function in x and k as in (40), 



d2 9 p -i,.6„4 da{£N 



£'ttN) 



(2tt 



d 2 q dx d?k 



47T 3 x(l — x) 



N,N 



[b;x,k) 



(53) 
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Thus the impact parameter distribution may be considered for fully exclusive (as well as fully inclusive) cross-sections. 

In the case of the fi — > ^7 example considered in Section III C the impact parameter amplitude A^l ? (6; x, k) 

is given by (39). The corresponding expression (51) for S^ n (b;x,k) is most easily evaluated by directly Fourier 

transforming I^A^i'^ 2 (q; x, k) | 2 , 
1+2+1 1 

S^(b;x,k) = (54) 
_ , 2 J fc V 2 X (2)ru\ ■ |fc| cos(0 b - fc ) Kxjmb) | K (mb) - \mh K l {mb) \ f kb 



[(l-x) 2 m 2 +fc 2 ] 2 W (1 - :r) 2 m 2 + k 2 2tt(1 - a:) 4tt(1 - x) 2 J V 1- 

The three terms correspond, respectively, to the virtual photon interacting (i) with the initial muon in both A^ 1 and 
(A 11 ) , (ii) once with the intial and once with the final muon, and (Hi) twice with the final muon. The imaginary 
part can be seen to arise from the angular correlation between the lepton scattering plane (defined by b) and the 

relative transverse momentum k in the final state. Representative plots of bS^i^l Xl ^)/4e 2 a; are shown in Fig. 1. 



V. DISCUSSION 



The impact parameter analysis of virtual photon induced transition amplitudes and cross-sections appears to open 
a new window on hadron dynamics. It is complementary to parton distributions in longitudinal momentum, and more 
economical in using data at all g 2 , not being restricted to the leading twist (q 2 — > 00) contribution. The analysis 
can be applied to any final (and initial) state, allowing to study systematic dependencies on, e.g., the mass, relative 
momenta and flavor content of the state. The J + component of the electromagnetic current needs to be isolated for 
a simple Fock state picture. 

Only Fock states that are common to the initial and final states contribute to the transition amplitudes (17), which 
are determined by the overlap (15) of the corresponding wave functions. This interpretation requires [4, 5] a frame 
like (16) with q + — , where the photon does not create or destroy quark pairs. This is analogous to DIS, where a 
parton model interpretation is possible only in "infinite momentum" frames with q + < 0. 

The momentum pf = p + q of the final state depends on the photon momentum q. Relativistic invariance requires 
that the momenta of all hadrons in / be parametrized as in (49), with the relative momentum variables Xi, ki being 
independent of pf. It is possible to form superpositions of final states through weighted integrals over the Xi and fcj. 
In the case of two-particle (ttN) final states we may thus consider states of the form (24) with photon matrix elements 

(7rN( Pf ;^)\J + (0)\N(p))^ / / -— (x,k)(7r( Pl )N(p 2 )\J + (0)\N(p)) (55) 

The pion and nucleon momenta are defined by (20) and we may freely choose the hadronic wave function tyf(x,k). 
The Fourier transformed amplitudes (17) get contributions only from quarks at b q = b, with the initial nucleon 



11 



and final 7riV states localized at zero impact parameter. The Fourier transform of the squared amplitude (51) gives 
the distribution of the impact parameter difference between the photon interaction vertices in the amplitude and its 
complex conjugate. 

The transverse shape of the contributing Fock states reflects only the distribution of the quarks struck by the 
photon, not that of the other partons. For example, both compact valence (Brodsky-Lepage [16]) Fock states and 
non-compact (Feynman [17, 18]) states may contribute to the elastic form factors of the nucleon at large photon 
virtualities \q\. Both types of states will contribute at small b q , since the photon interacts only with the x — > 1 
quark of the Feynman states, whose impact parameter is close to the transverse center-of-momentum (fojv = 0) of the 
nucleon. 

We expect the impact parameter distribution in j*N —¥ nN to contract as a function of the relative transverse 
momentum k between the final pion and nucleon. Only compact initial nucleons should have an overlap with irN 
states with high fc, in analogy to the observed color transparency of high energy pions dissociating into exclusive jets 
with high relative momentum [19]. 

Large angle photoproduction cross-sections are consistent with constituent counting rules [20, 21] at surprisingly 
low energies. Thus a{^p —> n + n) [22] and o(pfp — > K + A) [23] are both found to be cx i^M at ®cu = 90° • Even 
cr(7-D — > pn) [24] and (7(7 3 He —¥ pp(n)) [25] obey the rules, scaling as E@ M . The simplest theoretical prediction is 
based on perturbative QCD, which requires that only transversally compact Fock states contribute at large angles. 
Data on the q 2 dependence of large angle electroproduction would allow to to measure the actual width of the impact 
parameter distribution. 

According to the present analysis all contributing Fock states are common to the initial and final states. This does 
not require that heavy quarks Q of final states such as KA and DA C are present in the Fock states of the initial 
nucleon. Annihilations like QQ q — > q imply that heavy quark final states have Fock components with only light 
quarks q. We would still expect the impact parameter distribution to contract with increasing quark mass, since both 
the creation and annihilation of heavy quarks has a short transverse range ~ 1 /rriQ ■ 

The present analysis is also applicable to high energy diffractive processes such as 7*iV — > pN. As may be seen from 
(28) the momentum fraction x of the p meson decreases with the CM energy W (W 2 — p 2 ) as x ~ (m 2 p + k 2 )/W 2 . 
The quarks in the Fock states of the p meson have x q < x, and it seems likely that the virtual photon mostly interacts 
with such low- a; quarks. To the extent that the diffractive amplitude has a dominantly imaginary phase and conserves 
helicity the Fourier analysis may be done at the level of the amplitude, giving the impact parameter distribution of 
quarks with small x. 
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Appendix A 



Here we give analytic expressions for the matrix elements of the J + current in the QED amplitudes £(i — > £' /K7 and 
£7^ — > £' fi + (i~ at lowest order in a. They serve to illustrate the general expressions discussed in the text. 



1. £fi I'fif 



We denote the photon matrix element in (37) by 



A & = ^ T (^(Pi,A 1 ) 7 (p 2 ,A 2 )|J + (0)| Ai (p,A)) 



(Al) 



where p + q = p\ +P2 = Pf an d we parametrize pi,P2 as in (20). The expressions corresponding to (38) for the various 
hclicity combinations in the frame (16) are, 



A^ h + + l{q- lX ,k) = 2 ev ^ 
A^l(q;x,k) = 2ex 



e ■ k 



e_ ■ (fc — (1 — x)q) 



(1 - x) 2 m 2 + k 2 (1 - x) 2 m 2 + (fe - (1 - x)q) 2 
e + ■ k e_|_ • (fe — (1 — x)q) 



(1 - x) 2 m 2 + k 2 (1 - x) 2 m 2 + (k - (1 - x)q) 2 



A^l + J(q\x,k) = V2em^(l-xf 



A*_T_\{q;x,k) = and A%$ a (q;x,k)= ^'"^(-g; -fe) 



1 1 

(1 - x) 2 m 2 + k 2 ~ (1 - x) 2 m 2 + (fe - (1 - x)q) 2 

</J7,-A 



(A2) 



where m is the muon mass and e\ ■ k = —\e lX ^ k |fe|/y2 with k = |fe|(cos sin The Fourier transform (17) gives 



A^H(b;x,k) = 2eVx 



e ■ fe 



A^ + J{b;x,k) = 2ex5 



(1 - x) 2 m 2 + fe 2 
e + ■ fe 



5 2 (b)- 



s 2 (b) 



2v / 2tt(1 - 



-Kx{mb) 



exp —1 



{l-x) 2 m 2 + k 2 2v / 2tt(1-x) 



-JCi(m6) 



exp — z 



fc • b 

1 - x 

k ■ b 

1 - x 



A^^}{b;x,k) = V2emy/x~ 



(1-x) 2 



1 



n ^2 2x /( b )-^»H 
(1 — x) 2 m 2 + k 27r 



exp —i 



k ■ b 

l-x 



A^ri(b;x,k) = and (b; x, fc) = ^ 7 A '--A 2 (b; a, -fe) 



^7,A 



1/17, -A 



(A3) 



where 6 = |b|. The amplitudes (A3) are given by the overlap (15) of LF Fock amplitudes, according to the analysis 
of (40)-(42). 

The amplitudes for /17 final states with a fixed impact parameter b^ of the final muon are (cf. (45)), 



A»XU* x > b '» 



4\/27r 2 



lem 3 

x 2 

4\/27r 2 



4^7. 



y/x(l-x) e'^ i^i(mb^) 
Vx~(l ~ x) 2 K (mb^ 



-1 + e iqb '» 



-1 + e lq - b '» 



A^iqix^) = and (g; X, b)j = 

The Fourier transform (17) in q then gives, for all helicities (cf. (46)) 

A^ 2 (b- X ^) = ^x^x)r x t^K) 



A^\\ 2 (-q;x,b'^ 



- S W(b) + 5W(b-b') 



(A4) 



(A5) 
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In accordance with (15) also the amplitudes (A5) are thus given by the LF wave functions ip^ X2 (x, b'^) of the muon 
in impact parameter space [15], which describe /i~(p + , = 0; A) — > [i~(xp + , b'^; Ai) 7((1 — x)p + , —b'^ (1 — x)/x; A2). 

The Fourier transformed i\i — > t'fjry cross section (53) is expressed in terms of the squared amplitudes (51), 
S^^ 2 (b;x,k) = J cPb q (b q ; x, k) A^\* (b q — b;x,k). Their analytical expressions are actually most simply 

found from the Fourier transform of [A^'^ (q; x, k)\ 2 , 



S^(b;x,k) 



S^(b;x,k) 



S^(b;x,k) 



Ae 2 x 



k 2 /2 



[(1 -x) 2 m 2 + k 2 } 2 



S^ 2 \b)- 



\k\ cos(0 6 - 



4tt(1 - x) 2 



K (mb) — ^mb Ki(mb) 



2n(l-x) (\- x fm 2 + k 2 
k ■ b 



Kx{mb) 



exp 



l-x 



x 2 S^(b;x,k) 



(i-xY 



2e 2 m 2 x 

^^(m^jexp^-i— 



[(1 - x) 2 m 2 + k 2 ) 2 
k ■ b 



5^ 2 \b) 



(l-x) 2 



7T (1 - x) 2 m 2 + k 



■K (mb) 



» (6; x, k) = and Sj££ (b; x, k) = S^--x 2 (b; x, k) 



(A6) 



2. l-y* ->• £'fj,-(j, + 



We denote the photon matrix element in (37) for £-f* — > I'fj, /i + by 

1 



■Km = ^<M-(Pi,Ai)M + (P2,A 2 )|J+(0)| 7 >,A)) 



(A7) 



where p + q = Pi + P2 — Pf- The initial photon has virtuality p 2 and helicity A — ±1. At lowest order, in the frame 
(16) with pi,p2 parametrized as in (20), 



— 2ex' 2 Vl — x 
l-x 



A^t\(q;x,k) = V2em^/x{l - x) 

+ 2 + 2 

A^t\(q;x,k) 

' 2 2 

A^+\(q;x,k) 

2^2 

^1'+i(q;x,fc) 



1 



1 



M 2 + (k - (1 - x)q) 2 M 2 + {k + xq) 2 
e + • (k — (1 — x)q) e + ■ (k + xq) 



M 2 + (fe — (1 — x)q) 2 M 2 + (k + xq) 2 
A^Z\(q;x,k) 

Jb T 2 2 

and A^ 2 (q;x,k) = A^~* X2 (-q;x,-k) 



(A8) 



where M 2 = m 2 — x(l — x)p 2 and e\ ■ k = —\e lX ^ k |fe|/y2. The interaction of the virtual photon with the fi is given 
by the first terms in (A8), and the interaction with the n + by the second term. The Fourier transform (17) gives 



A^:\(b- X ,k) 

A^t\(b;x,k) 

TO Q 



A^\{b;x,k) 



A/27T 

+ieM 
V2-K 
1 - x 



y/x{\ - x) 



M 6^ 



k ■ b 



ji^r cxp {^i-r 



An 



' A/ b> 



exp +1 



k- b 



A^'t\(b;x,k) 

2 2 



A^_\(b- Xl k) = and ^; A A 2 (6 ;a ;,fc)- ^ A ^ A2 (6; S ,-fc) 



fiftyX 



afj,fj,, — X 



x exp +z 

x z \ x 



(A9) 
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where b = |b|. The amplitudes for fj, fi + final states with a fixed impact parameter b' of the final /i (cf. (45) are 



A 



Mb' 



exp +i<7 • o — exp — z 



1 - x 



^22 ^ 

A^illiq-^b'J 



+ieM 

4V27T 2 ' 



M b'„ 



1 — X 



exp I +iq ■ b' I — exp I — z 



zgfr' M 

1-a: 



1 — a; + a : 



: ^+i(q;x,b;) 

^+1 (q; x, b'j = and ^ (q; z, 6^) = (-g; s, ^ 



The Fourier transform (17) in q then gives 

A^ 2 (b;x,b^ = ^xJl^)^ X2 (x,b^ 



S^ib-b^-S^ib+^b'^ 



(A10) 



(All) 



In accordance with (15), the amplitudes (All) are given by the LF wave functions ifix \ (x, b'^) of the photon in 
impact parameter space, describing 7*(p + , b 7 = 0; A) —> /j,~(xp + , b^; Ai) — x)p + ', — b^ (1 — x)/a;; A2) : 



^22 ^ 



em ,Mb> 
17^ M ~x ] 



+ieM 
A^2n 2 



1 -x 



-i<Pb 



1 — x' 



^t\_(^K) = and K'xMK) = W-mUmW 



(A12) 



where M 2 = m 2 — x(l — x)p 2 . The wave functions were evaluated using the rules for LF wave functions in momentum 
space given in [7], Fourier transformed as in (31). 
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